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Abstract 

Voltage  loss  in  the  bipolar  plate  (BP)  is  induced  by  in-plane  current  in  the  BP,  which  arises  when  the  distributions  of  local  current  density  over 
the  surfaces  of  adjacent  cells  are  different.  We  show  that  potential  of  BP  satisfies  Poisson  equation  with  the  right  side  proportional  to  the  difference 
of  local  current  densities  on  both  sides  of  BP.  Solution  to  this  equation  is  obtained  for  BP  between  two  hydrogen  cells  with  the  single  straight 
channels  and  ideally  humidified  membranes.  The  general  relation  for  voltage  loss  in  the  BP  is  derived. 

©  2006  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Fuel  cell  modeling;  Voltage  loss;  Bipolar  plate 


1.  Introduction 

Bipolar  plates  (BPs)  separate  individual  cells  in  a  fuel  cell 
stack.  Each  BP  distributes  reactants  over  the  cell  surface  through 
the  system  of  channels,  collects  current  produced  by  individual 
cell  and  transports  this  current  from  one  cell  to  another.  BPs 
account  for  about  80%  of  the  stack  weight  and  45  %  of  its  cost  [  1  ] ; 
for  these  values  alone  it  is  clear  that  BP  is  one  of  the  key  stack 
components. 

In  this  work  we  will  focus  on  a  third  function  of  BPs  in  the  list 
above,  i.e.,  on  their  duty  to  transport  current  from  one  cell  to  an¬ 
other.  In  stack  modeling  it  is  usually  assumed  that  current  crosses 
BP  only  in  the  through-plane  direction  (see  e.g.  [2]).  This,  how¬ 
ever,  is  the  case  only  when  the  distributions  of  local  current  over 
the  surfaces  of  adjacent  cells  are  the  same.  However,  due  to  var¬ 
ious  factors  (different  stoichiometries  of  feed  molecules  in  the 
cells,  local  ‘‘spots”  of  contact  resistance,  CO2  bubbles  in  DMFC 
etc.),  these  distributions  are  usually  different.  Local  balance  of 
current  in  the  through-plane  direction  then  is  not  fulfilled  and 
significant  amount  of  current  produced  by  the  cells  flows  along 
BP  (in-plane  current).  This  current  results  in  voltage  loss  in  BP; 
calculation  of  this  loss  is  the  subject  of  this  work. 

It  is  worth  mentioning  that  local  voltage  drop  between  the 
two  bipolar  plates  is  equal  to  the  local  cell  voltage.  Therefore, 
in-plane  currents  in  the  BPs  change  the  distribution  of  the  anode 
and  the  cathode  polarization  voltages  over  the  cell  surface.  In 
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other  words,  the  problem  of  voltage  distribution  in  the  BPs  is 
coupled  to  the  problem  of  overpotentials  distribution  in  the  cells. 
Rigorous  approach  to  stack  modeling  thus  requires  solution  of 
equation  for  BP  voltage  together  with  a  complete  problem  for 
individual  cells. 

Voltage  loss  in  the  BP  depends  on  the  distribution  of  plate 
potential:  the  more  non-uniform  this  distribution  the  larger  the 
voltage  loss  (see  below).  Generally,  to  calculate  voltage  loss  we 
have  to  calculate  plate  potential  first. 

Direct  approach  to  this  problem  is  as  follows.  Potential  of  BP 
obeys  three-dimensional  (3D)  Laplace  equation;  the  boundary 
conditions  for  this  equation  provide  the  distributions  of  local 
current  density  over  the  surfaces  of  adjacent  cells  and  the  absence 
of  current  through  the  face  planes  of  BP  (see  e.g.  [3],  where  2D 
analog  of  this  approach  is  used  to  calculate  voltage  loss  in  the 
segmented  current  collector  of  a  single  cell). 

Below  we  will  show  that  potential  of  a  thin  BP  obeys  2D 
Poisson  equation  with  the  right  side  proportional  to  the  differ¬ 
ence  of  local  current  densities  in  the  adjacent  cells.  In  the  case 
of  BP  separating  two  cells  with  the  single  straight  channels  the 
problem  is  reduced  to  ID  Poisson  equation.  Solution  of  2D  or 
ID  equation  is  much  simpler  than  managing  3D  problem;  in  a 
ID  case  the  analytical  solution  can  be  derived. 

2.  Voltage  loss  in  the  bipolar  plate 

2.7.  General  assumptions 

The  thickness  of  BP  (^1  mm)  is  typically  2  orders  of  magni¬ 
tude  smaller  than  its  characteristic  in-plane  size  (^10  cm).  The 
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Fig.  1.  Sketch  of  the  2-cells  stack  with  the  large  aspect  ratio  L  w. 


variation  of  potential  across  the  plate  can  thus  be  neglected  and 
we  may  safely  assume  that  in-plane  current  in  the  BP  induces 
the  variation  of  plate  potential  only  along  the  BP  surface.  For 
simplicity  we  will  ignore  the  details  of  BP  geometry  (channel/rib 
structure)  and  consider  rectangular  BP  of  a  thickness  hv. 

We  will  assume  that  in  the  membrane-electrode  assembly 
(MEA)  current  flows  only  in  the  through-plane  direction.  Our 
simulations  [5]  show  that  the  fraction  of  current,  which  flows 
“along”  the  MEA  is  usually  indeed  negligible. 


2.2.  Quasi- ID  stack 

To  derive  an  equation  for  BP  potential  and  to  calculate  the 
respective  voltage  loss  consider  first  the  simplified  system:  the 
2-cells  stack  formed  by  the  two  identical  cells  A  and  B  with  the 
single  straight  channels  on  both  sides  (Fig.  1).  In  this  stack  the 
length  of  the  bipolar  plate  L  is  much  larger  than  its  width  w, 
L  w  (Fig.  1).  Since  L  w,  we  may  neglect  non-uniformity 
of  current  along  the  y-axis  and  consider  the  variation  of  plate 
potential  only  along  the  z-axis  (Fig.  1).  Physically,  each  cell  in 
our  stack  is  quite  analogous  to  that  considered  in  [4]. 

To  understand  the  origin  of  voltage  loss  consider  the  follow¬ 
ing  problem.  Let  local  current  density  ja  in  the  cell  A  linearly 
decreases  with  z  and  jb  in  the  cell  B  is  constant  (Fig.  2). 

As  in  real  stacks,  the  total  currents  in  the  cells  A  and  B  are  the 
same:  w  JQL  ja  d z  =  w  JQL  jb  d z  =  JLw,  where  J  is  the  mean 
current  density  in  the  stack.  However,  since  ja(z)  ^  jb(z)  there 
is  no  local  balance  of  the  through-plane  current.  In  other  words, 
local  current  generated  in  the  cell  A  cannot  reach  load  resistor 
moving  only  in  the  through-plane  direction  (Fig.  2).  Part  of  ja 
must  be  redistributed  along  z  to  match  the  shape  of  current  in 
the  cell  B.  This  redistribution  occurs  in  the  BP  and  induces  the 
variation  of  BP  potential  along  z. 

To  calculate  voltage  drop  along  the  BP  we  note  that  it  is  the 
the  difference  ja  —  jb,  which  has  to  be  redistributed  along  z. 
Elementary  current  (ja  —  j’b)w  dz  entering  the  plate  near  the 
point  a  must  be  transported  along  the  plate  to  the  point  a' 
to  provide  the  balance  of  local  currents  at  a  and  a'  (Fig.  2). 
The  respective  elementary  potential  drop  between  a  and  a'  is 
8Vaa'  =  (O'a  -  jB)wfc)(2(zo  -  z)r),  where  2(zo  -  z)r  is  the  to¬ 
tal  resistance  of  the  plate  between  a  and  a\r  is  the  linear  resis¬ 
tivity  of  the  plate  (Q  cm-1)  and  zo  =  L/2. 

Since  ja  —  jb  in  Fig.  2  is  an  odd  function  of  z  —  zo,  potential 
of  BP  at  z  =  L  is 


0  0.2  0.4  0.6  0.8  1 
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Fig.  2.  To  the  direct  calculation  of  voltage  drop  along  the  bipolar  plate.  Local 
current  in  the  cell  A  linearly  decreases  with  z;  in  the  cell  B  local  current  is 
constant. 


z)dz 


It  is  convenient  to  introduce  dimensionless  variables 


z 


where 


R  =  Lrwr = 


L 2 
hva 


V 

V  =  — 

RJ 


(1) 

(2) 

(3) 


is  the  plate  resistivity  (Q  cm2).  Here  hv  is  the  plate  thickness  and 
a  is  the  in-plane  conductivity  of  the  plate  material  (£2_1cm_1). 

Calculating  dimensionless  form  of  integral  (1)  with  j'A  =  2  — 
22,  jb  =  1  (Fig.  2)  we  find 

v(l)  =  \  (4) 


In  the  odd  case  (Fig.  2)  we  know  in  advance  that  (ja  —  j'b)w  dz 
has  to  be  transported  from  a  to  the  symmetrical  point  a' .  In 
the  general  case  of  arbitrary  ja  and  jb  these  simple  geometrical 
arguments  do  not  work  and  a  general  equation  for  plate  potential 
has  to  be  derived. 


2.3.  Equation  for  bipolar  plate  potential 

Consider  the  voltage  drop  8V  in  the  plate  on  the  interval  dz 
near  the  point  z.  Clearly, 

8V  =  Ir8z  (5) 

Here  /(z)  is  the  total  current,  which  flows  along  the  plate  at  z  and 
rdz  is  the  resistance  (£2)  of  the  domain  dz.  Current  /(z)  is  the 
sum  of  currents  ja  —  Jb  entering  the  BP  at  all  points  between  0 
and  z: 

I(z)  =  w  [  O'a  -  ,/b)  (6) 

Jo 
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Note  that  here  jA  and  jb  are  the  algebraic  values,  which  are 
positive  when  directed  along  the  normal  to  the  BP  surface  and 
negative  otherwise. 

Using  (6)  in  (5)  we  find 


dv 

3  z 


7b)  dz' 


Differentiating  (7)  with  respect  to  z  we  come  to 
d2V  dr 

— T  =  wrUA  -  7b)  +  I—  (8) 

ozz  oz 


We  see  that  potential  of  linear  bipolar  plate  obeys  Poisson  equa¬ 
tion.  If  r  is  constant,  Eq.  (8)  simplifies  to 

d2V 

=  wr(jA  -  jb)  (9) 

The  boundary  conditions  for  (8)  and  (9)  are 


V(0)  =  0, 


z=0 


(10) 


The  first  condition  establishes  the  reference  point  for  V;  the  sec¬ 
ond  one  follows  from  (7)  and  expresses  the  absence  of  normal 
current  through  the  end  face  of  BP. 

In  dimensionless  variables  Eq.  (9)  transforms  to 


d2V  . 

~^2  —  7  A  —  7B  (11) 


With  (1 1)  in  hand  we  can  calculate  the  distribution  of  poten¬ 
tial  along  the  plate  for  the  currents  in  Fig.  2.  Using  jA  =  2  —  2 z 
and  jb  =  1  in  (11)  and  integrating  we  get 

i Hz)  =  5?2  -  3 Z3 


At  2  =  1  we  get  V(l)  =  1/6,  which  coincides  with  (4). 


2.4.  Potential  of  the  two-dimensional  plate 


Eq.  (9)  can  be  generalized  to  the  case  of  2D  plate.  In  this  Sec¬ 
tion  we  will  derive  the  respective  equation  in  a  more  straight¬ 
forward  way. 

Consider  small  element  of  the  bipolar  plate  (Fig.  3).  Note 
that  in  this  Section  the  BP  surface  is  located  in  the  x-y  plane. 
Suppose  that  the  local  currents  in  the  adjacent  cells  A  and  B  are 
directed  as  shown  in  Fig.  3.  Balance  of  currents  in  the  element 
reads  (Fig.  3) 

jEhp  Sy  -  jwhp  Sy  +  jNhp  Sx  -  jshp  Sx  =  (jB  -  7a)  Sx  Sy 

(12) 


Note,  that  positive  is  the  current  directed  along  the  normal  to  the 
surface  of  plate  element. 

Dividing  both  sides  of  (12)  by  5v3y  and  noting  that 
7 e  -  jw  =  Sjx,  jN  ~  js  =  Sjy,  where  jx  and  jy  are  x-  and  y- 
components  of  in-plane  current  density,  we  get 


djx  djy  _  7B  ~  7A 

dx  3  y  hp 


According  to  Ohm’s  law  jx  =  —odV/dx  and  jy 
for  the  plate  potential  we  thus  find 

dfv  92v  _  JA  -  JB 

dx2  3  y2  hpG 


—a  dV/dy ; 


(14) 


Since  wr  =  1  /(/zpcr),  Eq.  (14)  is  the  2D  generalization  of 
Eq.  (9). 

For  the  2D  plate  the  natural  space  scale  is  the  plate  thickness 
hv,  rather  than  w.  Introducing  dimensionless  variables 


Vcr 

Jhp 


we  finally  get1 


d2V  i  d2V  _ 

UJ  +  ~ 7A  _  JB 


(15) 


(16) 


Normal  current  through  the  end  faces  of  the  plate  is  zero;  we 
thus  have  the  following  boundary  conditions  for  Eq.  (16): 


Here  Lx  and  Ly  are  the  plate  sizes  along  the  x-  and  the  y-axis, 

/V 

respectively.  These  conditions  determine  V  within  an  additive 
constant;  the  latter  can  be  defined  assigning  arbitrary  potential 

/V 

to  any  point  of  the  plate,  e.g.  taking  U(0,  0)  =  0. 


2.5.  Voltage  loss  due  to  in-plane  current  in  the  bipolar  plate 


To  derive  the  general  expression  for  voltage  loss  in  the  BP 
we  replace  the  BP  with  the  equivalent  resistor,  which  dissipates 
the  same  electric  power  and  transports  total  current  JS  in  the 
system,  where  S  is  the  BP  surface  area. 

Consider  first  the  quasi- ID  plate;  total  electric  power  dissi¬ 
pated  in  the  plate  is 


(17) 


Voltage  loss  on  the  equivalent  resistor  is  thus 


Uoss  — 


1 


*L 


rJwL  J o 


3V 
3  z 


d  z 


Note  that  j  =  j. 


03) 


1 
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In  dimensionless  variables  (2)  the  latter  equation  transforms  to 


(18) 


In  view  of  (7),  which  in  dimensionless  variables  is  dV/dz  = 
JqO'a  ~  Jb)  dz',  Eq.  (18)  is  equivalent  to 

fioss  =  £  (JTOa  - 


]'B)dz  dv 


(19) 


Thus,  in  a  ID  case  the  voltage  loss  can  be  calculated  not  solving 
the  problem  for  BP  potential.  Note,  however,  that  numerical 
solution  of  Eq.  (11)  and  subsequent  calculation  of  integral  (18) 
can  be  faster  than  the  calculation  of  double  integral  (19). 

In  2D  case  Eq.  (18)  takes  a  form 


Eloss 


dS 


(20) 


2.6.  Voltage  of  bipolar  plate  between  two  ideally  humidified 
hydrogen  cells 


To  illustrate  the  capabilities  of  this  approach  consider  local 
current  jx(z)  in  the  single-channel  hydrogen  cell  with  the  ideally 
humidified  membrane  [6,7] 


Mz )  =  h 


where 

A  =  -A'„(l-I) 


(21) 


(22) 


and  A  is  the  oxygen  stoichiometry.  It  is  easy  to  verify  that  jy 
dz  =  1. 


0  0.2  0.4  0.6  0.8  1 


z  /  L 


Fig.  4.  Local  current  densities  in  two  hydrogen  cells  with  ideally  humidified 
membranes  running  under  oxygen  stoichiometries  =  1.3  and  A-b  =  3.  To 
reach  load  current  Jb  —  Ja  must  be  transported  from  zone  a  to  zone  a'  through 
the  bipolar  plate  (thick  grey  arrow). 


Let  the  cells  A  and  B  in  our  quasi- ID  stack  operate  at  A  a  and 
Ab,  respectively  (Fig.  4).  Let  the  inlets  of  all  channels  be  located 
at  z  =  0  and  the  outlets  at  z  =  L.  The  local  currents  in  these  cells 
then  are  ja  =  jxA  and  jb  =  jxB .  Note  that  this  problem  cannot 
be  solved  using  simple  geometrical  arguments,  since  ja  —  jb  is 
not  an  odd  function  (Fig.  4). 

Using  jxA  and  j\B  (21)  in  Eq.  (11)  and  solving  we  find 

1 

Aa 


v(X)  = 


Aa 


ln(l  -  (1/Aa)) 


1  - 


Ab 

ln(l  -  (1/Ab)) 


+  C^B 


^a)z 


(23) 


Voltage  loss  due  to  in-plane  current  in  the  plate  is  given  by 
Eq.  (18).  The  difference  A  a  —  Ab  in  the  stack  is  usually  small; 
let  Ab  =  A  and  Aa  =  A  +  8X.  Calculating  integral  (18)  with  V 
(23),  substituting  Ab  =  A  and  Aa  =  A  +  8X  into  the  result  and 
expanding  over  8X  we  find 


Vi 


OSS 


1  + 


6A(1  -  X)  +  (1  -  2A)/ln(l  -  \/kj 


4(A  - 1)2  A2 


(<U)2  (24) 


where  f\  is  given  by  (22). 

The  function  in  the  square  brackets  in  (24)  is  shown  in  Fig. 
5.  Under  fixed  8X  the  voltage  loss  dramatically  increases  when 
A  tends  to  1.  Note  that  depending  on  the  feeding  scheme  the 
cells  in  a  stack  may  operate  at  different  stoichiometries.  Typi¬ 
cally,  the  cells  located  far  from  the  oxygen  inlet  operate  at  lower 
stoichiometries  than  those  at  the  inlet.  Eq.  (24)  shows  that  the 
voltage  loss  in  the  BP  between  these  “remote”  cells  may  be  sig¬ 
nificant.  Two  cells  running  at  different  stoichiometries  close  to  1 
form  a  “bottleneck”  for  current  transport  through  the  stack;  this 
bottleneck  may  reduce  the  overall  stack  efficiency.  Note  also 
that  the  voltage  loss  is  proportional  to  the  square  of  8X ;  it  is  thus 
important  to  keep  a  variation  of  A  between  two  adjacent  cells  as 
small  as  possible. 

With  A  =  1.1  and  8X  =  0.2  Eq.  (24)  gives  Vjoss  —  1.7  x 
10-2.  The  conductivity  of  carbon  BPs  o  varies  between 
102^_1cm_1  [8]  and  103£2-1cm-1  [9]  .  Taking  o  = 


Fig.  5.  The  function  F x  in  square  brackets  in  (24). 
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102  £2_1  cm”1 ,  L  —  10  cm,  hv  =  0.1  cm  and  7  =  1  A  cm-2,  we 
get  Vioss  —  140  mV.  This  value  is  quite  significant  for  fuel  cell 
applications. 

It  should  be  emphasized  that  Eq.  (21)  is  valid  when  the  polar¬ 
ization  voltage  of  the  cathode  side  ij  is  constant  along  z.  Voltage 
drop  along  the  bipolar  plate  may  violate  this  condition.  A  rig¬ 
orous  approach  requires  to  consider  the  problems  in  the  bipolar 
plate  and  in  the  MEA  simultaneously.  The  results  of  this  work 
will  be  published  elsewhere. 

3.  Conclusions 

We  have  shown  that  in-plane  current  in  the  bipolar  plate  arises 
when  the  distributions  of  local  current  over  the  surfaces  of  ad¬ 
jacent  cells  are  different.  This  in-plane  current  results  in  a  vari¬ 
ation  of  potential  V  along  the  BP  surface.  V  satisfies  Poisson 
equation;  the  right  side  of  this  equation  is  proportional  to  the 
difference  of  local  current  densities  in  the  adjacent  cells.  Ana¬ 
lytical  solution  to  the  problem  is  obtained  for  the  BP  separating 
two  ideally  humidified  hydrogen  cells  running  at  different  oxy¬ 
gen  stoichiometries.  This  solution  shows  that  the  voltage  loss  in 
the  BP  dramatically  increases  when  oxygen  stoichiometries  in 


both  cells  are  different  and  close  to  1 .  The  general  relation  for 
the  voltage  loss  in  the  bipolar  plate  is  derived. 
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